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Abstract. We review the status of the proton charge radius puzzle. Emphasis is given 
to the various experiments initiated to resolve the conflict between the muonic hydrogen 
results and the results from scattering and regular hydrogen spectroscopy. 


1 The proton charge radius puzzle 

The historical route to the proton charge radius (r p ) is from elastic electron-proton scattering. In 
a completely complementary fashion, it has been obtained also from “high-precision” laser spec¬ 
troscopy of hydrogen (H). Since a few years, “high-sensitivity” laser spectroscopy of muonic hydro¬ 
gen (/rp) offers a third way. The value extracted from /rp with a relative accuracy of 5 x 10 4 is an order 
of magnitude more accurate than obtained from the other methods. Yet the value is 4% smaller than 
derived from electron-proton scattering and H spectroscopy with a disagreement at the Ter level EH3- 
In the last five years as summarized in jg, |7[| various cross checks and refinements of bound- 
state QED calculations needed for the extraction of r p from /rp have been performed, together with 
investigations of the proton structure. Several suggestions in the field of “beyond standard model” 
BSM physics have been articulated, re-analysis of scattering data have been carried out and new 
experiments have been initiated. Despite this, presently the discrepancy still persists and the resolution 
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of the proton radius puzzle remains unknown. In this article, we summarize mainly the ongoing 
experimental activities which hold the potential to unravel the proton radius puzzle. 

As the atomic energy level are slightly modified by the nuclear finite size, it is possible to deduce 
the nuclear charge radius by performing spectroscopy of the atomic energy levels. In leading order, 
the energy shift caused by the nuclear finite size is 

A/ffinite size = ^I0 2 (O )|^ (1) 

where 0(0) is the wavefunction at the origin in coordinate space, m r the reduced mass, a the fine 
structure constant and n the principal quantum number. Re is the charge radius of the nucleus defined 
in a covariant way as the slope of the electric form factor ( Ge) at zero momentum exchange Q 1 


Re = -6 


dG E I 

dQ2\Qi = Q- 


( 2 ) 


Non-relativistically, Re is the second moment of the electric charge distribution pe of the nucleus 
R 2 e ~J dr p E (r)r 2 . 

The mj dependence of Eq. ([T} reveals the advantages related with muonic atoms. As the muon 
mass is 200 times larger than the electron mass, the muonic wavefunction strongly overlaps with the 
nucleus ensuing a large shift of the energy levels due to the nuclear finite size. Thus, the muonic 
bound-states represent ideal systems for the precise determination of nuclear charge radii. 

The proton form factors can be obtained from unpolarized differential cross section measurements 
of electron scattering off protons. In the one-photon approximation the elastic cross section is 


dcr 

dQ 


dcr | 

dClIMott 


1 

1 + T 


{g 2 e (Q 2 ) 



(3) 


where the Mott cross section applies for point-like particles, G E and Gm are the electric and magnetic 
form factors of the proton, r = Q 2 /4M 1 and e _1 = 1 + 2(1 + r)tan 2 (9/2) are kinematical variables. 
Commonly, the Rosenbluth separation is applied to disentangle the charge from the magnetic con¬ 
tributions by using the angle-dependence at fixed Q 2 . Therefore, by measuring the differential cross 
section at various Q 1 and angles 6, one obtains Ge(Q 2 ) and Gm(Q 2 ) and via Eq. ({2} the radius. 


2 Muonic hydrogen and possible new physics explanations 

The CREMA collaboration has measured two transition in pp: from the triplet (2.S' |,^ 1 - 2 P^J 2 ) [Q]] 
and the singlet (2 S - 2 P^p 1 ) 0 ] 2S-states yielding a radius of r p = 0.84087(39) fm. More specif¬ 
ically the two measured energy splittings, from the triplet liv, and from the singlet liv s states, can be 
combined to obtain both the 2S Lamb shift Ei = AE(2S - 2P\/2) and the hyperfine splitting £hfs: 


E l = ^hv s +^hv,- 8.8123(2) meV, (4) 

Ehfs = hv s -hv, + 3.2480(2) meV, (5) 

where the numerical values follow from reliable proton-independent corrections of the 2P states. Ex¬ 
perimentally n 

£“ p = 202.3706(23) meV = 48932.99(55) GHz (6) 

limited by statistics while the systematics effects is at the 300 MHz level. From theory fl] 

£* = 206.0336(15) - 5.2275(10)r 2 + 0.0332(20) meV (7) 
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where r p is given in fm. The first term accounts for QED contributions, the second for finite size 
effects, and the third for the two-photon exchange (TPE) contribution. From the Lamb shift an im¬ 
proved r p value free from uncertainties related with the HFS splitting has been determined, and from 
the HFS the Zemach radius, albeit not with the same accuracy as the charge radius. 

The consistency of the two /rp measurements, represents an important cross check of the muonic 
results. The typical systematic effects affecting the atomic energy levels are substantially suppressed 
in /rp due to the stronger binding. The internal fields and the level separation of the muonic atoms are 
greatly enhanced compared to regular atoms making them insensitive to external fields (AC and DC 
Stark, Zeeman, black-body radiation and pressure shifts). Thus /rp turns out to be very sensitive to the 
proton charge radius 0/z] 5 -dependence) and insensitive to systematics which typically scales as ~ 1 /in,-. 
The possible involvement of weakly bound three-body systems in the muonic hydrogen spectroscopy 
experiment @] has been ruled out by three-body calculations and by the experimental non¬ 
observation of sizable additional line broadening, line splitting and event rate decrease. 

The bound-state QED corrections which give rise to Eq. © have been computed by several 
groups OH 01 as summarized in and updated recently in Il3[ frill . Particular attention has been 
devoted to the TPE contribution which has been computed in two frameworks: one making use of 
dispersion relations and measured inelastic structure functions of the proton [15], the second from 
chiral perturbation theories E2]. Both ways provide consistent predictions. In the dispersion-based 
approach a subtraction term in form of an integral from zero to infinite Q 2 is necessary which can not 
be fully constrained by data. At intermediate Q 2 indeed a modeling of the proton is necessary Gl. 
The large majority of the community, agrees with a value for this subtraction term which is two or¬ 
der of magnitude smaller than the measured discrepancy of 0.3 meV H8-2(J. Even if improbable, 
in principle still a very un-smooth and physically unmotivated proton structure could be constructed 
which could shift the //p transition to explain the measured discrepancy 1210 . However, the published 
subtraction functions proposed to solve the proton radius puzzle would affect through the Cottingham 
formula the proton and neutron masses by 600 MeV 1220 which is quite implausible when compared 
with the measured and computed neutron-proton mass difference of 1.29 MeV 12311 . 

Explaining the discrepancy by this subtraction term can be interpreted as an exotic hadronic effect. 
Other effects occurring in the vicinity of the nucleus, as the breakdown of the perturbative approach 


in the electron-proton interaction 1241251. the interaction with sea fj ' // and e + e pairs [26] etc., have 
been suggested and need to be further investigated but at moment are not yet conclusive. Several BSM 
extensions have been proposed but the vast majority of them have difficulties to resolve the measured 
discrepancy without conflicting with other low energy constraints. Still some BSM theories can be 
formulated but they require fine-tuning (e.g. cancellation between axial and vector components), 
targeted coupling (e.g. preferentially to the muon or to muon proton) and are problematic to be 
merged in a gauge invariant way into the standard model 127 - 29h. 


3 Hydrogen experiments 

In a simplified way, the hydrogen S-state energy levels can be described by 

E(nS) = % + ^f (8) 

n A n J 

where R = 3.289 841 960 355(19) x 10 15 Hz is the Rydberg constant and approximately 


L ls - 8171.636(4)+ 1.5645/4 MHz 


(9) 
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the Lamb shift of the ground state ( r p expressed in fm) given by bound-state QED calculations. The 
different n-dependence of the two terms in Eq. © permits the determination of both R m and r p , from 
at least two transition frequencies in H, assuming Eq. <(9}. 

Being the most precisely measured transition ( u r = 4x10 15 ) 13C| 31 ] and because it is showing the 
largest sensitivity to the Lamb shift contributions usually the 1S-2S transition is used. By combining 
it with a second transition measurement, it is possible to determine the proton charge radius. When 
taken individually, the various r p values extracted from H spectroscopy by combining two frequency 
measurements (2S-4S, 2S-12D, 2S-6S, 2S-6D, 2S-8S, 1S-3S as “second” transition) are statistically 
compatible with the value from /rp. Only the value extracted by pairing the 1S-2S and the 2S-8D 
transitions is showing a 3cr deviation while all the others differ only by < 1.5cr. 

So the 4.4<x discrepancy between the proton charge radius from /rp and H spectroscopy emerges 
only after an averaging process (mean square adjustments of all measured transitions) of the various 
“individual” determinations and consequently is less startling than it looks at first glance. A small 
systematic effect common to the H measurements could be sufficient to explain the deviation between 
/rp and H results. This fact becomes even more evident if we consider the frequency shifts (absolute 
and normalized to the linewidth) necessary to match the r p values from /rp and H as summarized for 
selected transitions in Table Q] Obviously the discrepancy can not be solved by shifting the 1S-2S 


Table 1 . Transition frequency shift needed to match the r p values from H and /rp, expressed also relative to the 
stated experimental accuracy cr, and to the transition effective linewidths r eff . 


Transition 

Shift rel. to uncertainty 

Absolute shift 

Shift rel. to effective linewidth 

//p(2S-2P) 

100 cr 

75 GHz 

4T eff 

H(1S-2S) 

4'000 o- 

40 kHz 

40T eff 

H(2S-4P) 

1.5 cr 

9 kHz 

7 x io~ 4 r e ff 

H(2S-2P) 

1.5 cr 

5 kHz 

7 x io -4 r e ff 

H(2S-8D) 

3 cr 

20 kHz 

2 x io- 2 r eff 

H(2S-12D) 

1 cr 

8 kHz 

5xio- 3 r eff 

H(1S-3S) 

1 cr 

13 kHz 

5 x io -3 r e ff 


and the /rp measurements because it would require displacements corresponding to 4000cr and lOOcr, 
respectively. Expressing the required frequency shift relative to the linewidth as in the last column 
allows to better recognize some aspects of the experimental challenges. For example a shift of only 
7 x 1CT 4 T of the 2 S - 4 P transition would be sufficient to explain the discrepancy. A control of the 
systematics which could distort and shift the line shape on this level of accuracy is far from being 
a trivial task. Well investigated are the large line broadenings owing to inhomogeneous light shifts 
which results in profiles with effective experimental widths much larger than the natural linewidths |3|] . 

Another exemplary correction relevant in this context, named quantum interference, has been 
brought recently back to attention |[32h . and has lead to various reevaluations of precision experi¬ 
ments f33il . An atomic transition can be shifted by the presence of a neighboring line, and this energy 
shift 6E, as a rule of thumb, amounts to 113211 y 1 ~ - where D is the energy difference between the two 
resonances and T the transition linewidth. Thus, if a transition frequency is aimed with an absolute 
accuracy of T/x, then the influence of the neighboring lines with D < xT has to be considered. The 
precise evaluation of these quantum interference effects are challenging because they require solving 
numerous differential equations describing the amplitude of the total excitation and detection pro¬ 
cesses from initial to final state distributions and because it depends on experimental parameters such 
as the angular position and polarization sensitivity of the detectors, the laser intensity, direction and 
polarization, the initial population distribution among states, etc. 
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Generally speaking, transition frequencies involving states with large n are more sensitive to sys¬ 
tematical effects caused by external fields. Emblematic is for example the n 7 -dependence of the DC 
Stark effect. Motivated by the possibility that minor effects in H could be responsible for the observed 
discrepancy, various activities have been initiated in this field: 


• 25 — 4Pi/2 and 2S — 4 .P 3/2 at MPQ Garching & They are aiming at improving previous mea¬ 
surements by a factor of 5 down to an accuracy of few kHz which would yield a r p with less than 
2% accuracy when paired with the 1S-2S transition. Preparation of the 2S state by means of opti¬ 
cal excitation and an almost 47 t Lyman-alpha detection system are key elements to control the line 
pulling due to quantum interference on the 1 x 10 4 E level of accuracy required. 


15-35 transition at LKB and MPQ [35, 36]: The 2010 results from the LKB group delivered 


the second most precise transition frequency measurement in H with a total uncertainty of 13 kHz 
corresponding to a relative accuracy of 4.5 x 10 l2 . The error budget was dominated by statistics 
(12 kHz) and uncertainties in the velocity distribution of the atomic beam (3 kHz). A 1% accuracy 
of the proton radius will require a measurement of the 1S-3S transition with accuracy of about 
2 kHz. To reach this goal, the Paris group is presently pursuing the measurement of the 1S-3S 
transition at 205 nm wavelength using cw spectroscopy along the same line of investigations as 
in previous experiment. Special emphasis is devoted to the velocity dependent systematic effects. 
Oppositely, the MPQ group, to circumvent the difficulties related with the generation of the 205 nm 
light, has devised an experiment which uses pico-second frequency comb pulses. 


• 25 — 2P classical Lamb shift in Toronto Ml : The measurement of the 2S-2P energy splitting alone 
can lead to r p . Indeed, as this transition frequency does not depend on K m there is no need to 
combine it with a second transition frequency measurement. Microwave spectroscopy based on the 
Ramsey method of separated oscillatory field is used for this purpose. A factor of 5 improvement is 
anticipated, which implies a determination of r p to the 0.6% level. To reach this ambitious goal the 
position of the line has to be determined with 1 part in 10 4 . 


The “second” (beside the 1S-2S transition) transition frequency measurement in H can be inter¬ 
preted as a measurement of the Rydberg constant. An alternative way to an independent determination 
of the Rydberg constant, is to perform optical spectroscopy of H-like ions between circular Rydberg 
states where the nuclear size corrections are basically absent, the QED contributions small, and the 
linewidths narrow B38I1 . or via spectroscopy of positronium and muonium [39]. 


4 Scattering experiments 

The Mainz A1 collaboration at MAMI has measured in 2010 1422 precise relative e-p cross sections 
in the low-Q 2 regime (0.0038 GeV 2 to 0.98 GeV 2 ) and a wide range of beam energy and scattering 
angles |^1]. Two spectrometers were moved with overlapping angle settings while a third spectrometer 
was kept fixed and used as a luminosity monitor. 

As data are available only down to a minimal Q 2 , extrapolation to Q 2 = 0 is required to deter¬ 
mine the proton charge radius. Numerous works have been concerned with the issues related with 
this extrapolation procedure. The Mainz group found a satisfactory goodness of fit (x 2 = 1.14 for 
1422 points) through the use of flexible fitting functions (splines and polynomials) and the resulting 
radius reads 0 r p = 0.879(5) s tat(4) sys t(2) mo d e i(4) gl - OU p fm in agreement with the CODATA06 value of 
0.8768(69) fm based mostly on atomic measurements. It is important that the fit function be flexible 
enough to adequately reproduce the data, without being so flexible that over-fitting occurs. One so¬ 
lution to this problem is given in 00 where the low-Q 2 behavior of the form factor is constrained 
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by using large-/' assumption of the charge distribution. Reanalysis of the world data using these con¬ 
straints yields r p = 0.879(11) fm |42]. Another approach to address issues of over- or under-fitting 
data is the use of bounded polynomial z expansion (after conformal mapping) and constraining the 
expansion coefficients to decrease “perturbatively” with increasing order [40]. 

Conformal fits to the form factors were performed in [430 yielding r p = 0.870(23)(12) fm in 
agreement with the CODATA value. Another group fitting the 1422 Mainz data points found r p = 
0.840(15) fm [44] in agreement with the PSI result but with a^ 2 = 1.4 (using more flexible functions 
the same group found x 2 = 1.1). Noteworthy, the value of r p ~ 0.84 fm from analysis of scattering 
data using dispersion relations and vector-mesons dominance models was obtained by the same group 
prior to the publication of the muonic result. The use of a form factor model (for all Q 2 -ran ge) based 
on dispersion relation and vector-meson dominance introduce rigidity in the model which results in 
the larger^ 2 . So tension exists between the use of a physically motivated model giving a poorer fit or 
very flexible fit functions without physical constraints yielding better^ 2 . This tension could arise by 
inappropriateness of the theoretical model, insufficient treatment of experimental details, or due to an 
underestimation of the scattering cross section uncertainties. 

A wide-ranging study of possible systematic effects of the 2010 Mainz data has been recently 
reported in Special attention was devoted to the extrapolation procedure, to normalization fac¬ 
tors needed to smoothly combine the various spectrometer settings and to refinement of the radiative 
corrections. This reanalysis yields r p - 0.895(20) fm. When applied to the world data (excluding 
Mainz 2010) r p = 0.918(24) fm is found Even though some inconsistencies between data sets 
were found, which have led to the increased uncertainty of the r p extracted from scattering, it remains 
difficult to reconcile the scattering results with the muonic results. The only sure conclusion is that 
analysis of low-Q 2 scattering data is not simple and remains a matter of discussion. Because data at 
still lower Q 2 would be beneficial, two electron-proton experiments have been initiated: 


• PRad experiment at Hall B in JLAB & This experiment planned to operate at Q 2 down to 2 x 
1CT 4 GeV 2 aims to obtain r p with sub-percent accuracy. The experiment is based on a windowless 
target and a downstream calorimeter which allows to extend the cross sections measurements to 
smaller scattering angles. The need to measure relative cross sections at about 0.2% level requires 
knowledge of the angle to 10 /rrad accuracy which makes this experiment very challenging. 


• Initial state radiation at MAMI, A1 collaboration /’46f/: Making use of the initial state radiation 
techniques, where the initial electron momentum is degraded by photon emission, the momentum 
transfer to the proton is reduced. The scattered electron is measured with the usual spectrometers 
but no information on the photon is observed. However, by comparing measurements with Monte 
Carlo simulations accurate form factors can be determined down to Q 2 of 2 x 10 4 GeV 2 . 


Other scattering experiments can provide very important informations: 

• MUSE: muon scattering experiment at PSI & The MUSE experiment plans to measure p + - p 
and p~ - p as well as e + - p and e~ - p scattering down to 0.002 GeV 2 . By comparing negative with 
positive charges they will have an handle on the insidious TPE contribution. Comparison between 
electron and muon cross sections allows the elimination of common systematical effects, including 
some extrapolation uncertainties. Thus, this experiment has not only the potential to measure the 
proton charge radius absolute value to 2% accuracy, but also a possible difference between radii 
extracted from electron and muon scattering down to a relative accuracy of 1%. In this way possible 
muon-specific interactions can be disclosed. 

• Deuteron scattering at Mainz New e-d scattering data have been collected in Mainz, with the 
aims to extract a new value of the deuteron charge radius and break-up informations which may be 
used to compute the deuteron polarizability contribution in muonic deuterium. 
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Of relevance in this context is also the planned measurement at JLAB [49] aiming at the electric form 
factor of the mirror nuclei 3 He and 3 H to extract their charge radii difference, and the TREK program 
at J-PARC that scrutinizes K-decays to search for BSM physics motivated by the r p puzzle. 


5 Muonic deuterium and muonic helium ions 


In 2009 the CREMA collaboration measured two 2S-2P transitions in muonic deuterium (/yd). The 
2S-2P energy splitting in //d was determined with about 1 GHz accuracy, which corresponds to a 
relative accuracy of 20 ppm and 5% of the linewidth. 

Evaluation of the most challenging systematic effect, the quantum interference effect has been 
recently completed i33t l. Due to the proximity of two 2P states (4T apart), the quantum interference 
effects might be considerable. However, a quantitative evaluation of this effect, when accounting for 
the used excitation and detection schemes, the detector geometry, the laser direction and polarization 
etc., yields a line shift 6E < 0.001T, thus far below the statistical accuracy of our experiment. 

Moreover, the theory in muonic deuterium has only recently converged 52] to a state which 
allows a precise determination of the deuteron charge radius from the Lamb shift measured in /id. An 
impressive progress as been achieved in recent years both on the “purely” QED sector 112. 50l 5J, ], as 
well as in the computation of the TPE contribution 1 531-55], yielding [52]: 


AE pd (2S -2P m ) = 228.7766(10) - 6.1102(3) r 2 d + 1.7091(200) meV 


( 10 ) 


where the first term represent basically the “pure” bound-state QED contributions, the second term 
with r,i in fm the leading finite size contribution (including mixed radiative-finite-size corrections) 
and the third the TPE contribution. 

Combining our measurements with the prediction of Eq. ( ITOt we will obtain />/ with a relative 
accuracy of u r = 4x 10 4 limited by the TPE contribution. A second route to a precise r,y value is to 
combine the r p extracted from /up with the H-D isotopic measurements of the 1S-2S transition 0. A 
comparison of these two numbers will check the consistency of the muonic results and will give new 
constraints to BSM theories, e.g. if and how the "new force carrier" can couple to the neutron. 

In 2013 and 2014 we have measured for the first time 2 transitions frequencies in /r 4 He + and 
3 in yu 3 He + with relative accuracy of about 40 ppm. The uncertainty of the transition frequency 
measurement is entirely given by statistical uncertainty, since systematics effects or uncertainty related 
with the laser frequency calibration are < 2 ppm. To extract the nuclear charge radii from these 
measurements the corresponding theoretical predictions have to be known. Preliminary values read 

AE4 He (25 - 2P 1/2 ) = 1668.669(20) - 106.340 r 4 2 He + 9.52(30) meV (11) 

A£b Hc (2.S’ - 2P m ) = 1644.658(20) - 103.508 r 2 He + 14.66(40) meV . (12) 


Final numbers can not be given at this stage because it requires sorting out of the several contributions 
calculated by various groups using different frameworks [12, 50L 5jl 56]. Similar to fi d, in recent 
years there has been an impressive progression of the TPE predictions using state-of-the-art nuclear 
potentials, rendering the inelastic nuclear contribution with 5% accuracy Q5_7] . Such an accuracy opens 
the way to alpha-particle and helion radii determination with a relative accuracy better than 1 x 10 3 


which will be compared with the very precise value available from scattering |58[. Still, an important 
contribution of the TFE in /r 4 He + and in /r 3 He + , related with the intrinsic nucleon polarizabilities, has 
not yet been addressed by the community. “Simple” scaling as used in //d 05511 probably does not 
apply in this situation because of the smaller separation between nuclear and nucleon energies. 

Besides providing insights into the r p puzzle these nuclear radii represent benchmarks to check 
few-nucleon ab-initio calculations, or vice versa to fix low-energy coefficients (e.g. cu or ce of the 
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three-nuclei-interaction) describing the nuclear potential in effective field theories [59]. Moreover 
they can be used as anchor point for the 6 He- 4 He and 8 He- 4 He isotopic shift measurements [60]. 
The radii extracted from //He 1 measurements will be used to disentangle the 4cr discrepancy between 
two 3 He- 4 He isotopic shifts measurements |62|. and their knowledge opens the way to enhanced 
bound-state QED tests for one- and two-electrons systems in “regular” He + 163 1 and He 


6 Conclusions 

Various attempts have been made to find a solution of the proton charge radius puzzle which has been 
exposed by the Lamb shift measurement in muonic hydrogen. A plethora of theoretical works has 
been devoted in refining and rechecking the underlying theory necessary to extract the charge radius 
from the muonic measurements, in proposing theories beyond the standard model, in reanalyzing 
scattering data, and investigating proton structure. After all, the puzzle still persist. As a next step, 
new experimental inputs are required to provide guidance. Understanding of nuclear effects will be of 
primary importance for the interpretation of the next muonic measurements. 
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